Excitations in confined helium 
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We design models for helium in matrices like aerogel, Vycor or Geltech from a manifestly mi- 
croscopic point of view. For that purpose, we calculate the dynamic structure function of 4 He on 
Si substrates and between two Si walls as a function of energy, momentum transfer, and the scat- 
tering angle. The angle-averaged results are in good agreement with the neutron scattering data; 
the remaining differences can be attributed to the simplified model used here for the complex pore 
structure of the materials. A focus of the present work is the detailed identification of coexisting 
layer modes and bulk-like excitations, and, in the case of thick films, ripplon excitations. Involving 
essentially two-dimensional motion of atoms, the layer modes are sensitive to the scattering angle. 



I. INTRODUCTION 

There is much current interest in understanding the 
properties of superfluid He in confinement. In partic- 
ular, collective excitations of superfluid helium confined 
to silica aerogel have been studied by neutron scatter- 
ing since the early 90's Q, @, and by now a wealth of 
information about helium in aerogel. Vycor and Geltech 

has been collected HSHHHHI!lll[I3- Aerogel 
is an open gel structure formed by silica strands (Si02). 
Typical pore sizes range from few A to few hundred A, 
without any characteristic pore size. Vycor is a porous 
glass, where pores form channels of about 70 A diameter. 
Geltech resembles aerogel, except that the nominal pore 
size is 25 A [ll|. When liquid helium is placed in such 
an environment, it will first be physisorbed to the free 
surfaces of the matrices. Such an adsorption occurs in 
layers, the first layer of 4 He is expected to be solid; on 
a more strongly binding substrate, such as graphite, one 
expects two solid layers. In this work we are restricted 
to the liquid state, and the layer of solid helium is con- 
sidered as part of the substrate. 

Recent very accurate instruments have, contrary to 
earlier findings, revealed that the energies and lifetimes 
of phonon-roton excitations for confined 4 He are nearly 
equal to their bulk superfluid 4 He values. Specifically, 
the roton linewidth was found by Anderson et al. [T3] 
to be less than the instrumental resolution, 0.1/im. As 
expected, differences usually appear at partial fillings. 
The appearance of ripplons is tied to the existence of a 
free liquid surface; neutron scattering experiments show 
clearly their presence in on adsorbed fil ms |l 3Ul4l and in 
aerogel containing few layers of helium [jjj [lg ■ 

Since nothing spectacular could be seen in the bulk- 
like mode, the focus shifted to nearly two-dimensional 
"layer modes". The existence of collective excitations 
below the roton minimum has been demonstrated first 
by Lauter and collaborators (fjj Ea, UjJ |2(| ; these modes 
were identified as longitudinal phonons propagating in 
the first liquid layer close to the substrate |2lL I22I |23| . 
The nature of these modes has been found to be very 
resilient and quite independent on the substrate strength 
[2~H| . only for very weakly attractive substrates like Cs, 



the "layer phonons" acquire a transverse component and 
become similar to ripplons. Layer modes in helium have 
been observed in Vycor (2j| and in both fully filled and 
partially filled aerogel above wave vector 1.7 A -1 [2(|. 
Therefore, layer modes are the only excitations observed 
that are characteristic for helium films or helium in aero- 
gel or Vycor. 

The effect of disorder on macroscopic properties of 
both 3 He and 4 He has been investigated in detail, but 
there are only few theoretical studies on excitations that 
apply to confined helium. Locally, the density profile of 
the liquid can be assumed to be that of a helium film 
[27) or helium filling a space between two planes psj . 
Impurity induced disorder has been studied using Path 
Integral Monte Carlo j2^. At the moment there are no 
close-to-reality models of the confining structures com- 
bined with a proper theory of superfluid helium. Obvious 
simplified model systems are liquid helium filling random 
voids or helium between parallel planes. The former is 
computationally very demanding. Experiments have cov- 
ered the full range of helium thicknesses, from a bare 
monolayer to a filled sample, and the results show that 
layer modes are visible against the background of bulk 
excitations if there are about four or five layers of he- 
lium; in those cases we are basically dealing with a thick 
He film on Si. However, the pores in aerogel and Geltech 
are not filling uniformly, which has lead to the notion of 
connected sup erfluidity and localized Bose-Einstein con- 
densation 1 1 1| . Thus a well-defined phonon-roton spec- 
trum may persist above the critical temperature. 



II. THEORY 

The method of choice for studying the strongly inter- 
acting quantum liquid 4 He is the Jastrow-Feenberg the- 
ory. The theory is microscopic in the sense that it begins 
with the best available representation of the many-body 
Hamiltonian, 
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where V(\ri — | ) is the He-He interaction [3(J]. Com- structure function, 
pared to the bulk liquid, the essential new property that 
must be dealt with is the breaking of translational in- 
variance and anisotropy of the system. Such symme- 
try breakings may be induced by the substrate poten- 
tial U su b(ji), but they may also occur spontaneously like 
in the formation of quantum liquid droplets [3ll l32| . 
The ground-state theory has been throughly discussed 
in Ref. |53 and the corresponding excited-state theory is 
covered in Ref. I23L therefore we only need to review the 
key points. 

The ground-state wave function 
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is written in Jastrow-Feenberg form; the n-body corre- 
lation functions u n (ri, . . . ,r„) are determined by func- 
tional minimization of the energy. The theory yields the 
ground-state energetics and structure of the liquid and 
provides the raw material for investigating excitations. 

To introduce excitations to the system one applies a 
small, time-dependent perturbation that momentarily 
drives the quantum liquid out of its ground state. This 
causes the n-body correlations in the wave function pre- 
sented in Eq. to acquire time dependence, hence the 
excited state has the form 
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with the excitation operator 

SU(t) = y^fai(r^f) + y^Ju 2 (r l ,r j ;t) 
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The time-dependent correlation functions 
5u n {vu. . . , r n ; t) are determined by an action prin- 
ciple [33 
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where U ex t(t) is the weak external potential driving the 
excitations. 

Linearizing the equations of motion for the wave func- 
tion = l^o) + <5|*(i)) an d calculating the time- 
dependent component of the transition density 5p(r, t) = 
^o\p{ r )\8^(t)) +c.c. allows us to calculate the density- 
density response function ^(r,r';o;) defined via 



8p(r,w) = / rf 3 r' X (r,r»[/ ext (r». 



(6) 



Once the response function is known, one can apply 
the fluctuation-dissipation theorem to find the dynamic 



To obtain the dynamic structure function measured by, 
e.g., neutron scattering, one has to project S(r, r',u) 
onto plane waves: 

S(k,uj) = J drdr'e lk < r - r ^S(r,r'-,Lu) . (8) 

Up to this points the formulas are valid for any geome- 
try. Since the systems with slit or slab geometry under 
consideration here are translationally invariant only in 
the (ir — y) plane, but not in the z direction, only the 
momentum transfer k|| parallel to the film or slab is a 
good quantum number. In this case , we have 
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The truncation of the sequence of fluctuating correla- 
tions Su n in Eq. defines the level of approximation in 
which we treat the excitations. The excitation spectrum 
can be quite well understood [13 US Hi! by retaining only 
the time-dependent one- and two-body terms in the ex- 
citation operator The two-body terms £1*2 (ri, r^; t) 
describe the interaction of two excitations. The sim- 
plest non-trivial implementation of the theory leads to a 
density-density response function of the form p3| 

X(r,r»= (10) 

St 

where the (j)^(r) are Feynman excitation functions, and 
G st (uj) = [h[uj - uj s + ie]6 st + S st (w)] _1 (11) 

the phonon propagator. The fluctuating pair correlations 
give rise to the self energy |23| . 
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Here, the summation is over the Feynman states m, n; 
they form a partly discrete, partly continuous set due 
to the inhomogeneity of the liquid. The expression for 

(s) 

the three-phonon coupling amplitudes Vmn can be found 
in Ref. .23,. This self energy renormalizes the Feynman 
"phonon" energies oj n , and adds a finite lifetime to states 
that can decay to two lower-energy modes. The result- 
ing density-density response function has the structure 
of a Brillouin-Wigner (BW) perturbation formula. The 
approximate form of the self energy given in Eq. I|12l) is 
also closely related to the one obtained using the theory 
of correlated basis functions (CBF) [34L |3f| . As a final 
refinement to the theory, we scale the Feynman energies 
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ui n appearing in the energy denominator of the self en- 
ergy given in Eq. 1|12|) such that the roton minimum of 
the spectrum used in the energy denominator of Eq. Ijl2(l 
agrees roughly with the roton minimum predicted by the 
calculated S(k, ui). This is just a computationally simple 
way of adding the self energy correction to the excitation 
energies in the self energy itself. We shall use this scaled 
CBF-BW approximation for the numerical parts of this 
paper. 

A direct characterization of each mode is obtained by 
computing the transition density 

6p(v;t) = (¥ |p(r)|tf¥(t))+c.c. (13) 

and the transition current 

*j(r;t) = (*o|j(r) + cc (14) 

where p(r) and j(r) are the familiar one-body density 
and current operators. The transition density shows the 
density change (arbitrary amplitude) and the transition 
current the flow pattern of atoms in the mode. 



III. RESULTS 

We have modeled the confined quantum liquid by 4 He 
between two planar substrates. The distance between 
the substrate planes in our calculations is 40 A which 
is not too far from the diameter of aerogel strands and 
an intermediate value between the pore or channel diam- 
eters in Vycor (d ~ 70 A) and Geltech (d ~ 25 A). We 
assume translational invariance in a plane parallel to the 
surface. This is, of course, not exactly true in the above 
mentioned materials. However, our assumption should 
only change the details of the excitations spectrum at 
long wave lengths that are comparable to the pore size, 
whereas rotons arise from short range correlations. 

We have computed the dynamic structure function as 
a function of the scattering angle for a thick film and a 
filled gap. The density profiles for these two situations 
are shown in Fig. The interaction of 4 He particles 
with the walls is described by the usual 3-9 potential ob- 
tained from averaging Lennard-Jones potentials over a 
half space. This is valid if the walls are smooth and, 
hence, J7 su b(r) is a function of one coordinate only. Our 
approximation is legitimate because the lateral structure 
of surface is smoothed out by the first solid layer of he- 
lium atoms. The 3-9 potential used in this work is derived 
from the silicon-helium interaction |37|: we have supple- 
mented this potential by a 4-10 potential due to averaging 
Lennard-Jones 6-12 potentials for the 4 He- 4 He interac- 
tion over a plane. Thus, our substrate potential has the 
form 



with a well-depth D = 128 K and a range C3 
2000 KA 3 . The potential U 4 (z) has the form 



U suh (z) = U 3 (z - zO) + U 4 {z) 
where U 3 {z) is the common 3-9 potential |37| 



Uz(z) 



4Cf 



27D 2 



1 



(15) 



(16) 



U 4 {z) — A-Kepio 2 



10 



(17) 



Here, e = 10.22 K and a = 2.556 A are the usual deBoer- 
Michcls parameters for the helium-helium interaction. 
The areal density of the solid monolayer was taken to 
be p\ = 0.07 A 2 and a thickness zq = 3.3 A. The picture 
is consistent values for hectorite gaps [38| , the areal den- 
sity is somewhat lower than the one of the second solid 
layer on graphite [39j . 

The helium film considered here is thick enough to sup- 
port layer rotons; ripplons will also appear for the film 
model which has a free surface. In the filled gap case 
the amount of helium between the Si planes was chosen 
to correspond the filling of an aerogel sample surrounded 
by bulk liquid in equilibrium. The chemical potential of 
the confined liquid and that of a reservoir of bulk liquid 
turned out to be equal (p w —7.2 K) at n = 0.85 A~ 2 . 
For comparison, a 40 A thick slice of bulk liquid at the 
equilibrium density 0.02185 A -3 would correspond to the 
coverage n — 0.87 A~ 2 . 

Some excitations modes are bulk-like, these have no 
significant angle dependence. In the bulk liquid the roton 
energy decreases and the roton minimum moves to higher 
momenta as the density increases. Traces of this effect 
are also visible in the 3D roton: the roton energy is above 
it's equilibrium bulk value because of the low density 
regions between the layers. 

Layer modes, on the other hand, are entirely different. 
In films, layer rotons propagate mainly in the highest 
density liquid layer closest to the substrate, hence their 
motion is well confined to two dimensions and their prop- 
agation direction is always parallel to the surface. Their 
energy has a parabolic minimum at a fixed wave vector 
k|| g in the plane parallel to the surface. Upon chang- 
ing the orientation of the surface also the in-plane com- 
ponent of the total wave vector transfer changes. As a 
result, the location of the layer roton minimum shifts. 

We will show our results for the dynamic structure 
function as gray-scale plots. In these maps of the dy- 
namic structure function, darker areas correspond to 
higher value of S(k, ui). To emphasize low-intensity re- 
gions we use the scale S(k, uj) 1 / 4 , the same gray scale will 
be used throughout this work. To facilitate plotting we 
have also introduced a 0.05 K Lorentzian broadening of 
the structure functions. We measure the scattering angle 
from the substrate plane, so 9 = corresponds to grazing 
angle and = 90 degrees is scattering perpendicular to 
the substrate plane. 



A. Thick He Alms on silicon 

Fig. |21 shows the dynamic structure function of a thick 
4 He film on Si for scattering angles 8 = and 8 = 90 de- 
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FIG. 1: The density profiles of the film and the filled gap used 
in this work. The coverages of the film is 0.3 A -2 , and the 
filled gap has 0.85 atoms per A 2 . 



grees. The former shows a ripplon mode with dispersion 
hu>(k) ~ fc 3 / 2 , and also two layer rotons below the bulk 
roton are clearly distinguishable. Perpendicular scat- 
tering shows a discrete set of low-energy dispersionless 
modes, but some strong resonances appear also above the 
continuum limit —fi = 7.7 K. In the perpendicular case 
the familiar phonon-roton form is still recognizable, but 
beside it there are secondary maxima, which are roughly 
evenly spaced in k^. These strands appear because the 
projection onto plane waves in Eq.Q gives contribu- 
tions to multiple perpendicular momenta k^, and the 
dynamic structure function shows features coming from 
the Fourier transform of the density profile. A detailed 
description of how the strands appear was recently given 
inRef.|U 

Randomly oriented surfaces give rise to an angular av- 
eraged response. Fig.|3shows such a theoretical dynamic 
structure function, where contributions from several scat- 
tering angles between and 90 degrees have been added. 
While the contributions from bulk-like rotons add up to 
a single curve, the ones from 2D modes do not. The low- 
k region shows faints steps at about hui — K, 1.9 K, 
4.4 K, 6.0 K etc., reminiscent of the states turning disper- 
sionless as the angle increases; hence the fan-like struc- 
tures. The strongest of these comes from the ripplons, 
and is probably the only one that can be in experiments. 
As mentioned earlier, the layer roton minimum shifts to 
higher total momentum transfers as the scattering angle 
increases. This has not been seen in neutron scatter- 
ing experiments, which show a clear parabolic dispersion 
for the layer modes. The parallel direction gives the the 
strongest signal in neutron scattering, and we can only 
assume that the contribution from other angles in the 
raw data is simply too small to give any information to 
be inverted to S(k,Lu). 

Fig. 01 shows a fully angle-averaged dynamic structure 
function at k = 1.8 A -1 , near the roton minimum. Angle 




I 1 1 1 . 1 1 

0,0 0.5 1,0 US 2.0 2,5 3.0 

FIG. 2: The dynamic structure function of a thick 4 He film on 
a Si surface (coverage n = 0.3 A -2 ) at grazing angle (upper 
panel) and perpendicular to the substrate (lower panel). 



averaging causes the layer mode peaks to become asym- 
metric, while it has no effect on the bulk-like modes. The 
asymmetry is clearly visible in the first-layer peak, which 
has a broad tail on the high-energy side. This broaden- 
ing is just due to summing up of the multiple layer roton 
parabola, which show as separate curves in Fig. [21 A di- 
rect identification of the modes is given by the transition 
density and current, shown for the thick film in Fig.|Slfor 
6 = 0. The scale on the right corresponds to the transi- 
tion current, depicted as a vector field, for the duration 
of one oscillatory cycle which takes about 48 x TiwK -1 
picoseconds. The lowest energy mode for grazing angle 
is a layer roton, which involves atoms in the layer closest 
to the substrate. The next mode is loosely confined to 
the second layer, but subsequent modes spread through- 
out the liquid as there are no more well-separated layers 
available. From the third mode on the excitations are 
three dimensional, as well as they can be in a 15 A thick 
film. 

The situation in the perpendicular scattering is com- 
pletely different. The lowest panel of Fig. shows the 
mode at fko = 1.9 K, which already has one node in the 
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FIG. 3: The dynamic structure function of a thick film on a 
Si surface (coverage n = 0.3 A -2 ). The plot shows the sum 
of contributions at angles between and 90 degrees with 5 
degree steps. The dashed lines show the experimental layer 
roton dispersion for aerogel (upper curve) and Vycor (lower 
curve) by Plantevin et al. 10] ; The solid curve shows the layer 
roton and ripplon data from Ref. EH 
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FIG. 4: The fully angle-averaged S(k, u) of a thick film at 
k = 1.8 A" 1 . The coverage is 0.3 A" 2 . 



transition density. Particle number conservation forbids 
a nodeless mode, so the first real excitation is one where 
the liquid oscillates with respect to the substrate. 

Fig. El can be looked upon as a pictorial view of quan- 
tum evaporation. Higher energy modes involve motion 
of atoms near the surface, and finally modes with energy 
above minus the chemical potential (now fi = —7.7 K) are 
energetic enough to kick out atoms to continuum states. 
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B. 4 He filling space between two silicon surfaces 

Fig.0shows the energy and pressure in the case, where 
helium fills the space between the substrate walls. The 
equilibrium density of the filled gap by itself would be 



FIG. 5: Transition density (thick line) and transition current 
(arrows; right scale) for the low-energy modes at k = 1.8 A -1 
for a thick film on Si substrate and scattering parallel to the 
surfaces, 9 = 0. The grayscale plot shows the density profile. 
Excitations energies are indicated in the figures, the lowest 
panel corresponds to lowest energy. 
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FIG. 6: Same as Fig. E]for 9 = 90 (specular scattering). 
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FIG. 7: The total energy (solid line, left scale) and pressure 
(dashed line, right scale) of a helium filled 40 A wide gap as a 
function of coverage. Liquid in the middle of the gap becomes 
unstable at about 0.68 A~ 2 indicating the formation of two 
separate thick films on both walls. 



at about p = 0.78 A 2 , but the liquid is in balance with 
an external bulk 4 He reservoir at p = 0.85 A~ 2 . Below 
0.68 A~ 2 the liquid is no longer translationally invari- 
ant parallel to the walls, and the instability in the liquid 
solution is an indication of capillary condensation. At 
low density liquid forms two films covering both surfaces. 
In the filled gap case, if the density of helium between 
the walls is further increased, one observes a sequence 
of layering transitions: the number of helium layers in- 
creases. This filling scenario has been discussed in detail 
in Ref. H|. 

The (partially) angle averaged S(k,u>) in Fig.|S]is qual- 
itatively different from the film case result in Fig. [3] only 
in the low-momentum region. There is now only one 
strong maximum corresponding to bulk phonons on top 
of a faint step structure. The steps are a remainder of 
the discrete states in the perpendicular or nearly perpen- 
dicular directions. The angle averaged S(k, u) near the 
roton wave vector is quite similar to the one in the thick 
film case, depicted in Fig. 0] 

The transition densities and currents are shown in 
Figs. El and ^3 for scattering angles and 90 degrees, 
respectively. At grazing angle the low-energy excitation 
reside in the layers near the substrate, and pairs of modes 
corresponding to in- and out-of-phase oscillations are 
degenerate. The modes with even number of oscillations 
have their counterparts in the case of two films on both 
walls, as one can see by comparing Figs. El and The 
lowest mode depicted in Fig. EH with odd number of os- 
cillations have liquid oscillating back and forth in the 
middle of the gap, and with decreasing density this oscil- 
lation is finally able to divide the system into two films. 
The point where this mode becomes soft is the spinodal 
instability of the filled case. 
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FIG. 8: The dynamic structure function of the filled gap. The 
plot shows the sum of contributions at angles between and 
90 degrees in 5 degree steps. The calculated layer roton shifts 
to higher k as a function of angle. The dashed lines show the 
experimental layer roton dispersion for aerogel (upper curve) 
and Vycor (lower curve) by Plantevin et al. (lOf : the solid 
curve shows results from Ref. Ell 



IV. CONCLUSIONS 

In order to model 4 He confined in materials like aero- 
gel, Vycor or Geltech, we have calculated the energet- 
ics, the density profile and the dynamical properties of 
thick 4 He film on Si substrate and between two substrate 
planes. The simulation of confinement in spherical pores, 
which is perhaps the model that is closest to reality, is 
also feasible with our techniques, but computationally 
less efficient. For the excitations near the roton wave 
vector, where the most interesting phenomena are occur- 
ring, the use of a simplified planar geometry should not 
make much difference because the roton wave length is 
much shorter than any lengthscale on which the lateral 
structure of the aerogel changes. The two-wall geometry 
has actually been used in specific heat measurements by 
Gasparini and collaborators HH OH > but the dis- 
tance between the silicon planes was a couple of magni- 
tudes larger than in our model. For a large wall distance 
the mid section of the liquid is essentially unaware of the 
substrate. Thus one would have a thick region of basic 
bulk helium, which would show up almost exclusively in 
the dynamic structure function. 

The results for thick helium films are relevant for large 
pores at low filling, when opposite walls are approxi- 
mately noninteracting and the curvature of the liquid 
can be ignored. We have not considered the possibility of 
unevenly filled pores or a distribution pore sizes, because 
the characteristics of porous systems are very material 
dependent, and our aim was to discuss phenomena that 
are universal in such materials. 

There are presently two sets of data for the layer ro- 
ton, those of Refs. Holand ofl4ll When comparing theory 
and experiment, two aspects are of interest: One is the 
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FIG. 9: Transition density (thick line) and transition current 
(arrows; right scale) in the low-energy modes at k — 1.8 A -1 
in the helium filled gap at n = 0.85 A -2 for grazing angle. 
The grayscale plot shows the density profile. The transition 
density has arbitrary scale. Excitations energies are indicated 
in the figures, lowest panel corresponds to lowest energy. Due 
to the symmetric density profile there is pairwise degeneracy 
of modes corresponding to in-phase and out-of-phase oscil- 
lation between the left and the right side; we show here only 
the latter ones. 



energy of the layer roton, and the other one the disper- 
sion curve. As can be seen from the measurements of 
Ref. andliU shown in Figs. [3] and [HJ as well as from 
other measurements [2^, [2(| , that the energy of the layer 
roton depends sensitively on the nature of the substrate. 
This is expected because different substrates will lead to 
different densities of the first liquid layers and, hence, to 
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different phonon-roton spectra. Indeed, the dependence 
is so sensitive that it might be possible to determine the 
strength of the substrate potential from the energy of the 
roton gap. 



The second aspect is the shape of the dispersion curve. 
The measurements of Ref. EJshow a significantly smaller 
curvature of the dispersion curve than those of the other 
experiments, but they agree better with our calculations. 
Keeping the extreme difficulty of extracting these data 
from neutron scattering spectra in mind, we tend towards 
the view that the experiments of Ref. |4l| provide a more 
accurate description. This is so for two reasons: First, 
the curvature of the roton in purely two-dimensional 4 He 
is smaller than that obtained in Ref. ITfl and second, as 
shown in Figs. E3aiid|Hl the angular averaging leads to a 
further broadening. 
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scattering), except that now there is no degeneracy and we 
show all low-energy modes. 
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